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Abstract 

We study the effect of b —* sss scalar /pseudoscalar operators in B — > Krf'\ 4>K* decays. In 
the minimal supersymmetric standard model (MSSM), such scalar/pseudoscalar operators can 
be induced by the penguin diagrams of neutral Higgs bosons. These operators can be Fierz- 
transformed into tensor operators, and the resultant tensor operators could affect the transverse 
polarization amplitudes in B — > <j)K* decays. A combined analysis of the decays B — > (j)K* and 
B — » Kn^'\ including b — * sss scalar/pseudoscalar operators and their Fierz-transformed tensor 
operators originated from the MSSM, is performed. Our study is based on the followings: (1) 
Assuming that weak annihilations in B — > (j)K* is negligible and the polarization puzzle is resolved 
by Fierz-transformed tensor operators, it results in too large coefficients of scalar/pseudoscalar 
operators, such that the resulting B — > Kn^ branching fractions are much larger than observations. 
(2) When we take the weak annihilations in B — > <f>K* into account, the polarization puzzle can 
be resolved. In this case, new physics effects are strongly suppressed and no more relevant to the 
enhancement of the transverse modes in B — > <pK* decays. 



1 



I. INTRODUCTION 



Recent experimental results for polarization fractions in B ' — > 0(lO2O)i^*(892) ' are 



f L (B° - 0^*°) 



f L (B~ - 



f ± (B-^<f>K*-) 



0.506 ±0.040 ±0.015 


BaBar [1] 


< 0.45 ±0.05 ±0.02 




Belle [2] 


0.57 ±0.10 ±0.05 




CDF [3] 


0.227 ±0.038 ±0.013 


BaBar [1] 


< 0.31^ ±0.02 




Belle [2] 


0.20 ±0.10 ±0.05 

V 




CDF [3] 


j 0.49 ±0.05 ±0.03 


BaBar [4] 


1 0.52 ±0.08 ±0.03 


Belle [2] 


J 0.21 ±0.05 ±0.02 


BaBar [4] 


1 0.19 ±0.08 ±0.02 


Belle [2] 



(1) 



IAJ 2 , with 



Here, the polarization fractions fx (A = L, ||, ±) are given by f\ 
polarization amplitudes = A , An and A± being longitudinal, parallel and perpendicular 
modes in the transversity basis, respectively. Experimental results show that J'l ~ 0.5 and 
f± ~ f\\- On the other hand, the power- counting estimate in the standard model (SM) tells 
that the longitudinal mode is dominant [5]. In the SM, the QCD factorization (QCDF) 
calculation yields [5] f L : /y : f ± = 1 - 0(l/m 2 b ) : 0(l/m 2 b ) : 0(l/m 2 b ). The experimental 
results largely deviate from the intuition in the SM. Similar discrepancies have been observed 
in penguin-dominated B ±,a — > p ±,0 K*° decays [6, 7]. These discrepancies are referred as the 
polarization puzzle/anomaly in B — > (where V denotes a vector meson) decays. 

Solutions to the puzzle have been discussed within or beyond the standard model [5, 8, 9]. 
The recipe of fine-tuning form factors is proposed in [10]. Effects of final-state interactions 
are discussed in [11, 12]. Sizable annihilation effects are considered in [13-15]. As discussed 
in [14], the magnitude of annihilation correction is of O [l/ml log 2 m b /A h ^ . Furthermore, 
the effect is destructive to longitudinal, and constructive to transverse modes. Thus we may 
resolve the polarization puzzles by introducing annihilation effects. We note that, however, 
the perturbative QCD (pQCD) yields fi > 0.75 even with annihilation effects [8]. 

The b — > sg (where g denotes a gluon) operator, which enhances the transverse compo- 



nents, was discussed in [16]. However, it was found [13, 17] that the contribution due to the 
operator mainly affects the longitudinal mode. 

As for the solutions of the puzzle, the effects of NP-induced tensor operators are discussed 
in [17] and right-handed currents S7 M (1 + 7 5 )6g7 M (l ± 7 5 )g are in [18-20]. Because the 
right-handed currents may decrease the magnitude of \A \ and increase \A±\, it can explain 
the ratio However, the resulting \A\\\ <C \A±\ [18] is in contrast with the data 

\A\i\ ~ \A ± \. 

New physics (NP) contributions to B — > <pK* decays due to b — > sss tensor operators, 
first mentioned in [13], are systematically discussed in [17], and later the idea is applied to 
B — > pK* by considering the 4-quark tensor operators related to the processes b — > sdd and 

suu [21]. In the helicity basis, 1 four-quark tensor operators have leading effects to H (or 

H ++ ), but sub-leading to Poo- The possibility of solving B — > <pK* polarization puzzle by 
using four-quark tensor operators is extensively studied in [14, 17] and further investigated 
in [22-26]. 

In [17] the general approach of resolving the polarization anomaly of B — > <pK* by using 
four-quark NP operators is studied. There are two types of NP operators which are relevant 
to solve the polarization anomaly. They are tensor operators with er M!/ (l ±75) <E> a^(l ±75) 

structure. The tensor operator cr M!/ (l + 75) <E> o^iX + 75) results in Poo : H : H ++ = 

0(l/m b ) : 0(1) : 0{l/m 2 b ), while a^l - 75) ® a^(l - 75) leads to P 00 : H L_ : H++ = 
0(l/m b ) : 0(l/m 2 b ) : (9(1). 

The decays B — > PP (where P denotes a pseudoscalar meson) are sensitive to 
scalar /pseudoscalar 4-quark operators whereas B — > W are sensitive to tensor operators. 
Furthermore, it is known that scalar /pseudoscalar and tensor operators are not independent; 
scalar /pseudoscalar operators can be Fierz-transformed into tensor operators and vice versa. 
Therefore, the combined analysis of scalar /pseudoscalar and tensor operators for B — > PP 
and B — > VV modes will give more severe constraints about NP scalar /pseudoscalar and 
tensor operators. 

In this paper we focus on the b — > sss decay processes. We consider the 
scalar /pseudoscalar operators induced by Higgs penguin diagrams of the MSSM neu- 

1 Amplitudes in the helicity basis and the transversity basis are related by A n — H o, A\\ = (H ++ + 
H—)/V2, A ± = -(H++ -H — )/\/2. 
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tral Higgs bosons (NHB) [26-30], while the tensor operators which are obtained from 
scalar /pseudoscalar operators by the Fierz transformation. In this NP scenario, such tensor 
operators, contribute to the transverse polarization in B — > <pK* decays, whereas original 
scalar /pseudoscalar operators can affect B — > Krj^ decays. We study the consistency of 
both modes to see the validity of the scenario. One should note that this NP effect is further 
suppressed by m q /m s in the b — > sqq channel (with q = u or d) as compared with b — > sss. 
Although the recent observations of the sizable transverse fraction in B ± '° — > p ±,0 K*° decays 
may hint at large annihilation effects, the present study for decays B — > <f)K* and B — > Krf'^ 
can offer more severe constraints on the NP. 

The organization of the present article is as follows: In Sec. II, we summarize the for- 
mulation of MSSM-NHB scalar /pseudoscalar operators and its contributions to B — > Krf'^ 
and B — > <pK* decays. In Sec. Ill, we numerically analyze the decays for B — > Krf'^ and 
B — > (f)K*. Sec. IV is devoted to summary and discussions. 



II. FORMULATION 



A. SM and NP operators 



In the SM the effective Hamiltonian relevant to B — > Krj^ and B — > <pK* decays is given 



by 



Gp 

71 



v qh v* s {cMo\^) + cMo q M) 



,q=u,c 
10 



- VtbV t * s ^2ci(n)Oi(n) + c 77 (/i)0 77 (X) + c 8g {fi)0 S g{p) 



i=3 



h.c, 



(2) 



where the operators 0i=i,...,io are four-quark operators. Oj^ and Os g are electromagnetic 
and chromomagnetic dipole operators, respectively, fi is the renormalization scale. V q b and 
Vqs (<1 — u i c i t) are elements of Cabibbo-Kobayashi-Maskawa (CKM) matrix. The b — > sss 
four-quark NP effective Hamiltonian is given by 

26 



K 



NP 
eff 



= -%(W t :) J2 c iiv)Oi{n) + h.c, 

* i=U 



(3) 
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where Oi and q (i — 11, ... , 26) are four-quark NP operators introduced in [17], and corre- 
sponding Wilson coefficients 2 , respectively. Explicit forms oi Oi (i — 11, ... , 26) are shown 
in the following: 

(i) right-handed current operators 

On = Sa-y^l + 7 5 )6 Q %7^(1 + 75)5,3, 

(ii) scalar /pseudoscalar operators 

Ol 5 = S a (l + 75)&a + l5)S/3, 

17 = S Q (1 - 7 5 )6 a S/?(l - 75)S/3, 

Oi 9 = s Q (l + 7 5 )6 a ^(1 - 7 5 )s/3, 

2 1 = S a (l - 7 5 )6 a S/3(l + 75)S/3, 

(hi) tensor/axial-tensor operators 3 

023 = «a<V(l + 75^ 5/3^(1 + 75)S/9, 
25 = Sa^(l - l5)b a Sf3(J^(l - ^)sp, 

Since P — > PP (P — > VV) decays are not sensitive to the factorized tensor 
(scalar/pseudoscalar) 6 — > sss operators, it is a good approximation to use Oi with 

1 = 1, . . . , 22 (i = 1, . . . , 14, 23, . . . , 26) for P -> PP (P -> W). Some of these NP operators 
are not independent, and can be related with each other by the Fierz transformation: 

019 — — |Oi4, O20 = — |Oi3, O21 = — 106, O22 = — 1^5, 

23 = -40is - 80 16 , 24 = -80 15 - 40 16 , (7) 
25 = -40i7 - 80 18 , 26 = -80i 7 - 40 18 . 

Due to the Fierz transformation, we can introduce modified Wilson coefficients q, which are 
defined by 

Ci = Ci - l -c 3 with = (5,22), (6,21), (13,20), (14,19), (8) 

2 In [17], CKM factors and Wilson coefficients arc not separated. 

3 a^ v {\ ± 75)(T A " y (l =F 75) type operators Oi (i — 27, ... , 30) in [17] are found to vanish. 



012 = S Q 7^(1 + 75)6/3 5/37^(1 + l5)s a , 

Ou = s a 7^(l + 75)^/3 S/37 M (1 - 75)s Q , 

Ol6 = S a (l + 75)6/3 S/3(l + 75)Sq, 

Ois = s a (l - 75)6/3 s/j(l - 7 5 )s Q , 

O20 = S Q (1 + 75)6/3 S/3(l - 7s)s a , 
O22 = s a (l - 75)6/3 S/3(l + 7 5 )s a , 

24 = S Q <V(1 + 75)^/3 5/30^(1 + 7 5 )s a , 

2 % = s a a^(l - 75)6/3 S/36T^(1 - 7 5 )s Q . 
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and 



f-cA 




^15^ 










Cl6 


= M 


Cl6 




C18 


= M 


cis 


C23 




C23 




C25 




C25 


\C24J 




W/ 




\c 2 J 







/ 1 0-4 -8\ 

1-8-4 

with M = 

T2 -\ 1 

Hi l J 

Thus we can replace the Wilson coefficients by the effective ones: 

{ Cl ,c 3 } i = l,.. .,4,7,. ..,12,15,. ..,18, .7 = 5,6,13,14, 
for B — > PP decays, and 

{q,c,} i = l,...,4,7,...,12, .7= 5,6,13,14,23,. ..,26, 
for B — > VV decays, so that the decay amplitudes can be simplified. 



B. B — > ifr/W Decay Amplitudes 

In the SM, P> — > decay amplitudes are given by [31] 



where 



\/2T' p 



p=u,c 



A k v ^ [ 5 vu®2 + 2«3 + \al EW + 2/3g 3 ] 



4 2 Q; 3,£;U' 2 a i,Ew 



+ 03 + Ps,EW 



\/2T p 



+V2A Rv p [5 pu p 2 + a p 3 + 
+V2A Rr)(J) [5 pc a 2 + a p 3 ] 

+A V ^ K [ 6 vu{ui + fo) + «4 + 04,W +P3 + PIew] > 
+ V ^JO,<' ) K + "4 - \°%,EW ' \ a l,EW + Ps~ \Pl,EW + 

+V2A Rj) (,) [5 pc a 2 + af] 

[^4 ~~ l a 4,EW + 03~ \03,Ew\ ■ 

For the B — > PP decays ai t 2,3,4,3EW,4EW are defined as 

P P P P Pi Mo P 

"1,2 = Ol,2, «3 = a 3 - °5> «4 = °4 + r x X; 



a 3,£;VK — 



an — a 



P n P — n P I r M 2 P 
7) a 4,EW ~ u 10 ^ ' x 8' 
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with r% = 2m 2 K /m b (m q + m s ) and r^' = h 3 ^/ ' (f^ (!) m h m s ). A Mi m 2 are given by 

A MlM2 = • mlF^HO)/^. (16) 

Contributions from annihilation diagram are represented by (5q (Q = 3, 4, 3EW, 4EW, S3), 
where j3 Q = b Q ■ B MlM2 /A MlM2 with 

Gp Gp 

(where r = q or s) respectively. b 3 ^ 3 pw,4E\v are the coefficients due to weak annihilation of 
penguin operators. 653 is originated from the singlet penguin contribution which is intro- 
duced in [31, 32]. Note that following the approximation adopted in [31], we have neglected 
single weak annihilations /3si, Ps2, Ps3,ew, an d only keep ft $3- 

In the above results, we adopt aF i=1 ^ w and &q, given by the QCD factorization (QCDF) 
calculation in [31]. The NP effects due to scalar /pseudoscalar operators can be included 
by replacing c 5 and c 6 with the effective ones = c 5 ( 6 ) + Ac 5 ( 6 ) in the B — > Kr\ s decay 
amplitudes in the following way: 

Ac 5 = |(-ci 6 + ci 8 + c 20 - c 22 ), Ac 6 = |(-ci 5 + C17 + Ci 9 - c 2 i), for ol\ , (31 , (18) 

and 

Ac 5 = |(c 20 - c 22 ), Ac 6 = |(c 19 - c 2 i), for af , /3£, /3 53 . (19) 

Here C15, Ci6, cn, and cis are defined in (9). i] and 7/ mesons states can be regarded as 
mixed states of \i] q ) = -^(|-u-u) + |<ici)) and |?7 S ) = \ss) with a mixing angle (fi v [32]: 




cos^ -sin0 T; \ / |?7 9 ) 
sin <?!>,, cos(f) v J \\rj s 

Decay constants f q £ , pseudoscalar densities are defined by 

(?7 (/) (p)k"7 M 75?|0> = -^f Q v o)P^ (^ (,) (p)|s7 m 7 5 s|0) = -i/^p", 
2m g (r ? W(p)|g 75 g|0) = --^ w , 2m a <77W(p)|g 75 g|0) = 



(20) 



(21) 



with m q = (m u + rrid)/2. As for explicit forms of f^, h 9 ^ and form factors F ( f^' ?< ' (g 2 ), we 
summarize in Appendix A. 



C. B — > </>if* Decay Amplitudes 



Decay amplitudes of B — > can be decomposed as 

A(S - 07?*) = E ~Hfihi (22) 

where 

# = E W ( T ^,a + Vk\b) , (* = 0, ±), (23) 

P=M,C 

is the amplitudes in the helicity basis. Amplitudes for the emission topology (the T4 part) 
is given by 4 

E W^&M = (- V * V ts) { Ah (BR*,t)- [«S + °4 + 4 - 44 -\{a h 7 - r% + a\ + a h 10 )] 

Jl ~ A (BK*,4>)+ [°11 + fl 12 + a 13 ~~ r x°14] 

+^(Bi?*,<?>)T+ [ a 23 + | a 24] + ^(W*,c/>)T- [ a 25 + 2 fl 2e] } ' (27) 



p=u,c 



with given by 



r* = 2^!M. (28) 



The coefficient "1/2" in front of 2 6 can be realized as follows. Take O23 as an example. Because, under 
the Fierz transform, 2 3 can be written by 

O23 - (1/2)^(7^(1 + 75)6^ Bpa^il + J 5 )s a - 6* Q (1 + 75)^/3 + 75)s Q , (24) 

therefore, in the factorization limit, we obtain 

(^K*\0 23 \B) = (<P\sa^(l+ l5 )s\0)(K*\sa^(l+ l5 )b\B) 

+l(^\s a^(l+^)s a \O)(K*\s a a^(l + l5 )bp\B) 

= ( 1 + ^f) (01^(1+ l5 )s\0)(K*\sa^(l + l5 )b\B}. (25) 

Note that the second term in the right hand side of (24) gives no contribution since the local scalar current 
cannot couple to <f>. Similarly, in the factorization limit we have 



{4>K*\0 2i \B) = [ -L + ± ){cj>\sa^(l+ l5 )s\0)(K*\sa^(l + l5 )b\B). (26) 



The same procedure can be applied to the matrix elements containing O25 and 2 &- 
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Coefficients (i — 3, . . . , 10) have been calculated in QCDF [15, 31]. However instead of 
C5 and cq, C5 and c 6 should be used in the calculation of 05 and a& (see (11)). 

, Cj±l 



+ 0(a s ), with i = 23, 24, 25, 26, 



(29) 



where the radiative corrections are negligible. We summarized the explicit form of a\ for 
i — 11, . . . , 14 due to the right-handed four-quark operators in Appendix C. 
In (27), coefficients Af BVuVa)± and A h {BVi y 2 ^ T± are given by 



G 



(BK*,<(>)=f = ^(0(g^i(/^))|^(l-75) S |O)(^(p / ,e 2 (/ i ))|s 7 ,(lT75 



)b\B(p)) 



G F 



{if<f>m<p 



-2i 



m B + m K * 



e^*e v *p^V{q 2 ) 



(m B + m K *){e\ ■ e^A^q 2 ) - (e* • p)(e* 2 ■ p) 



2A 2 {q 2 



A h - - 



(BK*,4>),T± - ^2 



V2 

Tif<t>m<p 

%(0(g,e 1 (/ i ))|^ S (l±75)|O)(^(p / ,e 2 (^))|^(l±7 5 )fe| J B(p)) 



m B + ttik* 



h (30) 



^fli^.eTefp^Uq 2 ) 
T^T 2 (q 2 ) [(el-e* 2 )(m] 



m 



±8tT 3 (q 2 )(e* 1 -p)(e* 2 -p) 



K* 
2 

mi 



2 2 

771 £ TYlj^* 



2(el ■ p)(e* 2 ■ p)} 
h 



or in the explicit forms 



A ± 

(BK* ,4>) — 

A ± 

^{BK*,4,)+ 



G 

T-^iiUm^irriB + m x *)Mi(mJ) - 6A 2 (mJ)], 

(m B + m K *)M\mA =F — — — — V{rn 4 



V2 
Gf 



(if*™*) 



/ \ j / 9 \ BPc ... v 

-(to b + m^)Ai(mJ =F ■ V{m 4 



m B + m^* 



m B + 



4^^ )T± = T^m^rn 2 B [h2T 2 (ml) - h 3 T 3 (ml)}, 
A (Bif*)T + = --^4(</J')m|[±/ 1 r 1 (m;)-/ 2 r 2 (m;)], 
4^)t- = -- v |4(i/J)m B [±/ 1 T 1 (my + / 2 r 2 (mJ)], 



with a = (m B — — m 2 K *) / {T.m^rriK*) , b = (2m 2 B p 2 c ) /[m^mK*{jn B + m#»)] and 



(31) 



(32) 



/1 = 2p c /m B , f 2 = (m 2 B - m 2 K *)/m B 

,2 ^2 ^2 V^2 
1 w / ' ' ' ' ' ' 

h 2 = 



1 



{m 2 B -m\-m\*){m\-m 2 K *) 
4 p c 



4p>| 



Hi 



2mK*m < f ) \m B — m K * 



(33) 
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Here we have used decay constants and form factors denned in Appendix B. 

Weak annihilation contributions (the 7^-part) to B° — > (j)K*° and B~ — > <fiK*~ decay- 
amplitudes in helicity basis can be given by 

p=u,c 

E W^S?-^ = ^ {(" [6& + *W + b h 5 ] + k^c • b h 2 } , 



where 



and 



p=u,c 



h h 



N? L 



.G F 

B<t>K* = i~^=fBfK*f<t>, 



c^Af + c b {Af + A{' h ) + N c c 6 A f / 1 
c 9 A^ h + c 7 {Af + A(' h ) + N c c 8 A f 3 ' h 



b\ 



C F 



c u Af + ^(Af + 4' h ) + N c c l4 Al" 



with h — 0, — , +. Building blocks A l ^ ,h can be found in Appendix of [15]. 



D. Scalar and Pseudoscalar operators in the MSSM 



(34) 



(35) 



(36) 



In the MSSM, scalar and pseudoscalar operators can be induced by neutral-Higgs boson 
(NHB) penguin diagrams. We refer such operator as the MSSM-NHB scalar /pseudoscalar 
operators. In [28], b — > sit (where I denotes a charged lepton) scalar /pseudoscalar op- 
erators Q 12 induced by the MSSM-NHB penguin diagrams are considered, b — > sqq 
scalar /pseudoscalar operators can be obtained by replacing Higgs-l-£ vertex with Higgs-g-g 
vertex [29]. They are 

015 = s(l + 75 )6 V] ^-q(l + 75 )g, Oi6 = Si(l + -y 5 )bj V —q 3 (l + 75)9*, 
17 = s(l - 75)6 V —g(l - 75)^, 0i8 = Si(l - 75)^ V — 9j(l - 7s)%, 
0i 9 = s(l + 75 )& V - 75)9, 020 = Si(l + 75)^ V ^(1 - 75)*, 



m b m b 
q g 

021 = S(l - l5 )bJ2 — 5(1 + 75)9, 022 = Si(l - 75)^ E — ^ + 75)ft, 



(37) 



m b 
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where q — u,d, s, c. 5 The Wilson coefficients C;(/x) of Oi with i — 15, . . . , 22 at fj, — m w are 
given by [28, 29]: 

2 2 

C 15 (m w ) = ^(C s + C P ), C 17 (m w ) = ^(C' s - C' P ), 

2 2 

C 19 (m w ) = ^(Cs - C P ), C 21 (m w ) = ^{C' s + C' P ), (38) 

C i = (i = 16,18,20,22), 

and four-quark tensor operators are not directly induced. Here 

n (>) _ 4 g 2 s m 2 b cos a 2 + r s sin 2 a m~ g , , d LL(RR) ,<ilr(lr*) 



3Af # 2 sin 2 # w m^ cos 2 (3 m b ' 

V = — — ~r"( r p + tan2 p)—fb( x )^ LlRm ^3 RiLR,> > 



(39) 



with r s = m 2 H0 /m 2 h0 , r p = m 2 A0 /m 2 z0 , x = m\jm 2 g and At = V tb V t *. g 2 and g s are gauge 
couplings for the weak and strong interactions, respectively m h o, m H o and m^o are masses 
of neutral Higgs bosons hP, H°, A , respectively, a is the neutral Higgs mixing angle and 
tan (3 is the ratio of the two Higgs vacuum expectation values, m g and m^ are the gluino 
mass and common squark mass, respectively. Factors ^ L ■> $23 R an d ^33 R are down- type 
left-light second-third, right-right second-third and left-right third generation squark mixing 
parameters, respectively. The loop function f' b {x) is defined as f b (x) = (x 2 /2)d 2 f b o(x) /dx 2 
and fbo(x) is defined in [29]. f b (x) is given by 

Because Oi (i = 15, 22) and Oi are related by 

— O t C O it (41) 
m b 

the Wilson coefficients for Oj(/i) with i — 15, . . . , 26 at fj, — m w are given by 

m> Tfi 
ci 5 (m w ) = — Cis = D(A - B)£, c 17 (m w ) = -^C 17 = D(A - B)£, 
m b m b 

Tfl TYl 

ci 9 (m w ) = —C 19 = D(A + B)£, c 21 (m w ) = —C 21 = D(A + B)£, (42) 
m b m b 

Ci = for i = 16, 18, 20, 22, 23, 24, 25, 26, 



5 Precisely speaking, in the two doublet Higgs model, couplings of the light neutral Higgs h° to up-type 
quarks are suppressed by tan (3 compared with the down-type quarks. Therefore we can neglect contribu- 
tions of up-type quarks. 
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where 



A = 



D = 




B = 



to 2 . 




Since the Wilson coefficients for (sb)(dd) scalar and pseudoscalar operators are suppressed 



by md/m s , we neglect contributions for operators with q = d and we consider only b — > sss 
type operators. In (43), we note that D is almost positive and real because f' b (x) < 0, 
Re At < and the imaginary part of Xt is negligibly small. 

There are three enhancement factors: tan/3, 1/to^o and l/m H o in Eq. (43). Therefore, 
if R is sizable and neutral Higgses, A and H°, are sufficiently light and tan (3 is large, 
then the Wilson coefficients for scalar /pseudoscalar operators can be large enough. In the 
present paper, for simplicity, we neglect the mixing between scalar /pseudoscalar and tensor 
operators through renormalizat ion-group equations (RGEs) since such a mixing effect is 
small. Thus we assume Cj(m&) oc Ci(mw) and we use same symbols A, B, D, £ and £' at 
/i ~ m b . 

Under the existence of the MSSM-NHB scalar /pseudoscalar operators (43), Ac 5i6 , given 
in (18), (19), in B — > Kr/ S decays are rewritten by 




(44) 



where = |fW| exp(i0 ( - / - ) ). As for the Wilson coefficients in B — > decays, we use the 
replacements (11). They are given by 




(45) 



C5 = C 5 , Cn = C12 = C13 = 

c 23 = ^ J D(A- J B)£, c 24 = - 
c 25 = ±D(A - B)?, c 26 = 



\D{A-B)i'. 
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FIG. 1: (B - A)/B for various tan/3 and m^. 
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T77 



0.01 
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As for 023-26) we parametrized the following coefficients that appear in (27). 



°23 + 2°24 
1 

°25 + 2° 26 



1 £-4 

8iV c B 
1 B - A 



i(6±4>) 



DB\£\e 



DB\i'\S s '^\ 



(46) 



8iV c 5 

where the ct s corrections are negligible. However, we still parametrize the strong phases 5 
and 8' here. Actually the strong phases consistent with zero in the fit. Here and below the 
helicity labels are omitted for 023-26 since these coefficients very weakly depend on their 
helicities. 

The factor (B — A)/B depends on the details of neutral Higgs sector. In FIG. 1, we 
have plotted the value of (B - A)/B for various m A and tan (3 in the MSSM. In the MSSM, 
{B - A) jB is always smaller than one and -0.1 <{B-A)/B< 0.3 for tan/3 > 1. When 
(B - A)/B ~ 0.2, the ratio ((023(25) + |a 24( 26))/Ac 6 | w 0.02(5 - A)/B w C(10" 2 ). 



III. NUMERICAL ANALYSIS 



A. Numerical Inputs 



We summarize input parameters in Table I. As for B — > K* vector and tensor form 
factors, we follow the light-cone sum-rule (LCSR) results [33], defined as 



F{q 2 ) = F(0) exp( Cl g7m| + c 2 q 2 /m 2 B ), 



(47) 
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TABLE I: Input Parameters 



Decay constants [5, 31] f n = 131 MeV f K = 160 MeV f B = 210 MeV 

/, = (1.07 ± 0.02) U fs = (1-34 ± 0.06)/* 

f K * = 218 MeV /£, = 175 MeV / = 221 MeV /J = 175 MeV 
B meson parameter [31] \ B = 200±g 50 MeV 
7] — r]' mixing angle [31] 0^ = 39.3° ± 1.0° 

CKM parameters [35] A = 0.807 ± 0.018 A = 0.2265 ± 0.0008 \V ub /V cb \ = 0.881^;^ 

3 ( 7 ) = (67.4+2|)° sin 20! (sin 2/3) = 0.688±g;g|| 
B meson lifetimes [36] rgo = 1.530 ps r B - = 1.638 ps 
B -> P form factors [31] ^^(0) = 0.28 F^ K (0) = 0.34 



for F = Ao, Ai, A2, V, Ti, T2, T3, where c\ and C2 are listed in Table III of [17]. We use 
m B ± ps m B o = 5.279 GeV. We parametrize \ B to be m B /\ B = J* 1 dy&Bi(y)/y, where 
$si(y) is one of the two £> meson light cone distribution amplitudes with y being the 
momentum fraction carried by the light spectator quark in the B meson. As for the renor- 
malization scale we use p = m b /2. In the calculation of the hard spectator and the weak 
annihilation, we adopt ///, = y/A h • p ~ 1 GeV (with A h = 0.5 GeV) corresponding to the 
hadronic scale. 

We also include hadronic uncertainty parameters defined in [15, 34] : X^f , X^f (M = 
K,rf'\ K*), and X% *. For simplicity, in the fit, we assume that X^ H = X^ H = X^^ ^ in 
B — > Krj^ decays, and parameterize them as [34] 

^iT^^ + ^exp^^log^), 0<p A , H <l. (48) 

In B — > (f>K* decays we have fixed X L = m B /A h and X H = \og{m B / K h ) (i.e., p B = Pl = 0), 
because in these decays the branching ratios are very insensitive to them. 

As for NP effects, we use DB\£\, DB\£'\, 0, 0', 5, 5' and (B - A)/B as the independent 
parameters. We have constrained weak and strong phases to be |0^| < n and \8^\ < ir/2, 
respectively. In the fit, for simplicity we consider the two scenarios: (i) NP-(A) for which 
f = and (ii) NP-(B) for which f = 0. 
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TABLE II: World averag 


;cs of observables for 5 — > 


Krf'^ are shown in 


the second column [36- 


46]. Upper limits are at ! 


90% CL. In the third and fourth columns we have shown best fit values 


for combined fit with corrections received from B - 


-> annihilation. 


Corresponding best fit 


parameters are shown in 


Table IV and in FIG. 2. 


Best fit values of B — ► </>.ff* are shown in 


Table III. 








Observable 


Experiment 


Combined Fit with </>if* ann. 






NP-(A) 


NP-(B) 


B(B + -► rfK+) x 10 6 


69.7±|;f 


68.9 ±2.2 


68.9 ±2.2 


B(B° -► i]'K°) x 10 6 


64.9 ± 3.5 


66.3 ±2.1 


66.3 ±2.1 


B(B + -► 7?^+) x 10 6 


2.2 ±0.3 


2.2 ±0.3 


2.2 ±0.3 


B(B° -► r/K°) x 10 6 


(< 1-9) 


1.5 ±0.3 


1.5 ±0.3 




0.031 ±0.021 


0.031 ± 0.021 


0.030 ±0.018 




0.09 ± 0.06 


0.02 ±0.01 


0.02 ±0.01 




-0.29 ±0.11 


-0.30 ±0.11 


-0.31 ±0.10 


A CP (B° - r/K ) 


(N.A.) 


-0.20 ±0.11 


-0.14 ±0.07 


-ricpSicojiBO^r/K ) 


0.61 ±0.07 


0.70 ± 0.03 


0.70 ±0.01 


- V CpS k0t} (B° - r/if ) 


(N.A.) 


0.89 ±0.17 


0.73 ±0.03 



B. Experimental Data 

In the fit for B — > Krj^ decays, we use 7 observables including 3 averaged branching 
fractions, 3 direct CP violations and the B° — > K°rj' indirect CP violation —r] C pS Kri '. Here 
S Kj] (,) is defined by 



_ 2Im(A / ) q A(B° -> K^f) 



with q/p ~ e -2 ^ 1 for B° d) and r/^p is the CP eigenvalue of \Kg L rf^). The value of —rj CP S Kri (,) 
should be close to sin20i in the SM. The experimental data for B — > Krf'^ are listed in 
Table II. 

For the B — > <f>K* decays we have 20 observables, which include B°~ — > tpK*°~ branching 
fractions (B), polarization fractions (fi, f±), and CP asymmetries (Aq P , A° CP) Aq P ), phases 
of polarized modes and phase differences (A0y, A0j_), where A^p = \A\\ 2 — 
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Ea I^a| 2 )/(Ea I^aI 2 + Ea I^a| 2 ), A" cp = (|A A | 2 - |A A | 2 )/(|A A | 2 + |A A | 2 ), 6 A = arg(A A Mo), 
A(f>x = | arg(A A /A • Aq/Ax) with A = 0, ||,±. The experimental data for the B — > 
decays are shown in Table III. 

C. Combined Fits 

If we ignore the annihilation effects in B — > decays, the resulting x 2 m in ^ 170 is 
too huge; i.e. we cannot have a reliable fitting result. This is the fact that if the B — > 
<pK* polarization anomaly was mainly due to the tensor operators induced by the Fierz 
transformation, then the NP effects would lead to too large B — > Kr)^ branching ratios as 
compared with the data. 

Once the B — > (f)K* annihilation effects are included, we can see that the x 2 m m * s drasti- 
cally small. In Table IV, we have summarized the best fit values of the x 2 min an d parameters. 7 
The resulting 8 and 5' are consistent with zero, which are also consistent with the fact that 
the a s -corrections to the tensor operators are negligible. The fitted results for B — > Krj^ 
and B — > cj)K* are collected in Tables II and III. 

The results obtained in [17], where the weak annihilation effects are not included, are 
given by 8 

|S 23 ,dy| = 4.36^ x lO" 3 , |5 25 , DY | = 5.38l°f x lO" 3 , (50) 
to be compared with our present upper bounds, 

|a 23 | < 7.1 x 10" 4 , |a 25 | < 6.1 x 10" 4 , (51) 

which are extracted from FIG. 2 and (46), and are much smaller than the values in (50). This 
indicates that the contributions of tensor operators induced from the scalar /pseudoscalar 
operators in the MSSM-Higgs are too small to explain the polarization puzzle, while the 
puzzle can be accommodated by the weak annihilations effects. 

6 In BaBar measurements [49], instead of Aq P , the asymmetries of f^ 1 and f^ 1 are denned, f^ 1 and f^ 1 
are the polarization fractions measured in B and B decays, respectively. 

7 The errors of parameters in Table IV are obtained from the error matrix (covariance matrix) at the global 
minimum of \ 2 - The error matrix is the inverse matrix of the curvature matrix of chi-square function 
with respect to its free parameters. The errors of best-fit values in Tables II and III are estimated from 
the same error matrices for each NP scenario. 

8 Here a 2 3 and a 2 5 are defined by a 2 3 = a 2 3 + and a 2 5 = a 25 + ^a 2 6, respectively [17] . We have factored 
out the CKM factor. 
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TABLE III: World averages and best fit values of observables for B° — > 4>K*° (upper) and B — > 
(lower) [1-4, 36, 37]. In the third and fourth column, best fit values for combined fit of 
B -► Krf'} and B -► <£K* for NP-(A) and NP-(B) scenarios, including the contributions from 



5 — ► annihilations, 


are shown. Corresponding 


best fit parameters 


are show in FIG. 2 and 


Table IV. Best fit ^ 


/alues 


for B — ► Krf'} are shown in 


Table II. 




Observable 




Experiment 


Combined Fit with cf)K* ann. 








NP-fAl 


NP-fBl 


StotXlO 6 




y.o ± u.o 


O A _i_ n c 

• / . t: _1_ U.U 


q q _i_ n 6 






1 n n _i_ i i 
1U.U ± 1.1 




9 4 + 6 


h 




n /ioi 4- n oqo 
u.4yi ± u.UoZ 


n aqi + n 095 


n aqo 4. n noA 

U.^tfO _l_ u.u^t: 






U.OU ± U.UD 


0.499 ± 0.028 


0.486 ± 0.025 


f± 




n o^o -L- n hqi 
U.ZDz ± U.Uol 


n 047 _i_ n ni 

V . / _1_ U.U-LO 


n oon _i_ n noi 

U.aiJi? _l_ U.UZ;± 






n on _i_ n r\K 
U.zU ± U.UD 


n 04Q _i_ n ni n 


n 041 _i_ n 092 


A tot 
Jir' p 




_n ni _i_ n n« 

U.Ul _1_ W . W V J 


_n no i n 

U.UO ZIZ U.UtJ 


n no + n m 

u.uu m u.u± 






— U.U1 ± U.Uo 


02 + 03 


n no + n nn 


j4/~r p 




n no _i_ n nv 
U.Uz ± U.U ( 


_n f)3 ± 03 


nni + n 01 






U.l / ± U.ll 


ri fie _|_ n 07 


n nn + n nn 






n 1 1 4- n i o 
— U.ll ± U.1Z 


0.02 ±0.03 


0.00 ±0.00 






n oo _i_ n oc; 
U.zz ± U.ZD 


0.05 ±0.07 


0.00 ±0.00 


011 




9 o 7 +0.14 
z -°'-0.13 


2.34 ±0.09 


2.53 ±0.02 






2.34 ±0.17 


2.33 ±0.09 


2.52 ±0.02 






2.36 ±0.14 


2.49 ± 0.06 


2.56 ±0.03 






2.58 ±0.17 


2.48 ± 0.06 


2.55 ±0.03 






0.10 ±0.14 


0.02 ±0.06 


0.00 ±0.00 






0.07 ±0.21 


0.03 ± 0.08 


-0.01 ±0.00 


A</> ± 




0.04 ±0.14 


0.03 ± 0.07 


0.00 ±0.00 






0.19 ±0.21 


0.04 ±0.09 


-0.01 ±0.00 
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FIG. 2: Contour plots for A X 2 = X 2 ~ X 2 min in Re(DB^) v.s. Im(DB^) [or Re(DB£') v.s. 
lm(DB£')] for the NP scenario-A [or NP scenario-B]. Allowed regions of A% 2 < 1, 1 < A% 2 < 4 
and 4 < A% 2 < 9 are shown by dark, medium-dark and light-gray regions, respectively, "x" 
symbol indicates the location of the global minimum, x 2 min . The origin corresponds to the SM. 
The circle at the origin indicates the allowed upper-limit from the B s — > data. 
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TABLE IV: Best fit parameters obtained in NP-(A) and NP-(B) with considering B — ► (f>K* anni- 
hilations. Numbers with (*) indicates that they reach the upper or lower bound in the parameter 
space. Best fit values of B — ► Krj^ and B — ► are shown in Tables II and III, respectively. 
Weak phased NP parameters (DB£, DBg) are plotted in FIG. 2. 

Combined fit with (f>K* annihilation effects 





Scenario (A) 


Scenario (B) 


X 2 min/d-0.f. 


9.8/17 


15.5/17 


(B - A) IB 


0.55 ±0.76 (*) 


0.43 ±0.19 (*) 


5 , 5' 


S = +30 ± 85° 


5' = -38 ± 70° 




1.00 ±0.33 (*) 


0.51 ±0.41 (*) 


<f> A [Kv {,) ] 


117 ± 17° 


53 ± 47° 


PH[Kr}W] 


0.11 ±0.67 (*) 


1.00 ±0.59 (*) 


c/>h[Kv {>) } 


69 ± 156° 


-167 ±25° 


Pa[<PK*] 


0.57 ±0.04 


0.55 ±0.02 


4> a [4>k*} 


-96 ± 3° 


-85 ±3° 



D. Consistency with SM and B s — ► 



The current upper-bound for the branching fraction of B s — > ji + ji [37, 47] at 90% CL is 



B(B S -> < 7.5 x 1(T 8 . 



(52) 



The branching fraction of b — > s££ due to operators Of^ (with i = 7,9,15,17,19,21) for 
£ = e, /x, r 9 is given by 



x { (1 - 4m 2 ) 



r W _ r W 4- r W - r (e) 
c 15 c 17 "T c 19 c 21 



cS? + 4? - eg - 4? + 2m(c^ - cH ^3) 



9 We have defined o\ as O^L = s(l — 75)6^(1 ± 75 )f, and 0>^ (i — 15,..., 21) with replacements: 
s(l ± 75 )s -» £(1 ± 75K and s 7 ^(l + 75 )s -» (1 + 75 )s. 
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where rh = m^/m Bs , cf^ are the Wilson coefficients of Of^ at /x = m fc , and we have used 

(0|*7 5 6|B.) = ^ 3 . (54) 

Here we note that if RGE effects are not large, cf 1 ~ (m^/m s )q and 

(cib - en + ciB - c 21 ) = 2DA(£ - C) = 2 (1 - (B - A)/B) DB{^ - £'), 

(55) 

(cis + ci7 - ci9 - C21) = -2DB(e + eO- 
Using -B s 's lifetime tb s = 1.437 ps, mass m# s = 5.366 GeV, decay constant fs s = 215 ± 
25 MeV, quark masses m& = 4.9 ±0.1 GeV, and m s = 145 ± 25 MeV, we obtain upper 
bounds of DB^ as 

\DB£ {,) \ < 9.2 x 10~ 4 . (56) 

In FIG. 2 we have shown the SM (DB£ = DB£' = 0) and the upper bound of the NP 
effect constrained by the B s — > decay as the origin and a small circle at the origin, 

respectively. In the figures we have also drawn the contours of Ax 2 = 1,4 and 9, where 
Ax 2 = X 2 ~ X 2 mm- The B s — > fi + fi~ allowed region shares partly the Ax 2 < 1 (ler) region 
in NP-(B), and is just outside of the Ax 2 < 4 (2a) region in the scenario NP-(A). Since in 
both cases the B s — > data and the SM are located within contours where x 2 /d.o.f. is 

sufficiently small, we can safely conclude that our two scenarios are consistent with the data 
for B s — > decay and with the SM. 

In [23], the authors discuss the scalar /pseudoscalar operators induced by R-parity vio- 
lating interactions in the supersymmetric standard models. Because they did not take into 
account the weak annihilation effects and possible constraints from B — > Krj^'\ large con- 
tributions due to tensor operators to explain the B — > <pK* polarization puzzle are required 
and therefore the estimated magnitudes of the effects of scalar /pseudoscalar operators are 
much larger than the upper bound of B s — > ji + ji~ . 



IV. SUMMARY 



We have studied the scalar /pseudoscalar operators, and tensor operators where the latter 
are obtained from scalar /pseudoscalar operators by the Fierz transformation, in B — > <pK* 
and B — > Krf'^ decays. We have considered the scalar /pseudoscalar operators induced by 
penguin diagrams of MSSM neutral Higgs bosons. 
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Without the weak annihilations in B — > <pK*, we cannot obtain any reasonable solution 
to explain both B <pK* and B — > if^W decays simultaneously in the NP region (—0.1 < 
(5 — A) fB < 1) of the MSSM induced by the neutral Higgs bosons. Taking into account 
weak annihilation effects in B — > we obtain best fit results in good agreement with the 
B — > 1^77 W and B — > data. From the fitted parameters we estimate the magnitudes 
of the contributions due to NP tensor operators. They are, however, much smaller than 
the results of [17], which are introduced to explain the B — > <f)K* polarization puzzle. 
The polarization puzzle can be explained mostly by weak annihilation effect, as pointed 
out in [13-15]. The contributions of NP operators are constrained mainly by the fit of 
B — > Kr}W data. While our results may allow non-vanishing NP effects, the data for the 
decays B — > Krf'\ B — > 4>K* are consistent with the B s data as well as the SM 

prediction. 

Finally, we remark on the recently observed large longitudinal polarization fraction fi 
in B — > 0_K|(143O) [1]. If tensor operators play an significant role in B — > VT (where 
T denotes a tensor meson) decays, fi may significantly deviate from unity. The current 
B — > 0^(1430) experiment seems to be consistent with our conclusion since in our anal- 
ysis the effect due to tensor operators is found to be very small. However, in the present 
study we cannot exclude the possibility that sizable NP effects contribute directly to ten- 
sor operators, instead of scalar /pseudiscalar operators, and, moreover, a cancelation may 
take place between weak annihilations and contributions due to NP tensor operators in the 
B — ► 0fQ(143O) decay. For the point of view of the new physics, B — > 0fQ(143O) may 
be sensitive to the B — > tensor form factor which can be further explored from the 
B -> fsT* (1430)7 decay. 
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APPENDIX A: DECAY CONSTANTS AND FORM FACTORS FOR 77 AND 7/ 
MESONS 



The \rj) and \rf) meson states are defined as the mixed states of \i] q ) and \rj s ), as stated 
in (20) [32]. In this section we summarize the notations in [32]. Decay constants are 
given by 



(Al) 



(A2) 



f?, = fq C °S <t>r,i fv = ~f* Sm <^> 

fl> = f q sin (f) v , f °, = f, cos (f) v , 
and in the same way, pseudoscalar densities h^f are defined as 

h q v = h q cos (f) v , h s v = -h s sin <f) v , 
h q v , = h q sin (f) v , h s n , = h s cos 4> v , 

where h qyS are defined by 

h q = f q (ml cos 2 4> v + m 2 , sin 2 <p v ) - V2f s (ntf f - m 2 ) sin <fi v cos <f) v , 
h s = / s (m 2 , cos 2 <f) v + m 2 sin 2 (j) v ) - — / 9 (m 2 ; , - m 2 ; ) sin 0^ cos </y 

B — > 77 ^ form factors are defined as 

and in the present paper we take Fi = F B ^ 7r (0) and F 2 = 0. is the decay constant of the 
pion. 



APPENDIX B: DECAY CONSTANTS AND FORM FACTORS IN B — > VV DE- 
CAYS 

We have used 

(<£(g,e*)mo) 
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/^ej , (Bl) 
" W^***^(? 2 ), (B2) 



m B +m K *)£* K * lJi A 1 {q ) - (e* K * ■ Pb)(pb + Pk* 



^m B + m K * 



-2im K * £ -^^-q» [A 3 (q 2 ) - A (q 2 )} , (B3) 



for current operators, and 



(K*(p',e 2 )\sa^(l + l5 )b\B(p)) = ie^e^p'^T^q 2 ) 

+{4,(m 2 B - m 2 K .) - (s;-p)(p + p'),}T 2 (q 2 ) 



(B4) 



{4-Pb) 



% — 2 — —(p + p)v 



T 3 (q 2 ), (B5) 



for tensor operators. In (B3) and (B5), A 3 (0) = A o (0), 71(0) = T 2 (0) and 
M q 2 ) = mB+mK * Mq 2 ) - mB ' mK ' A 2 {q 2 ). 



Ivtik* 



(B6) 



APPENDIX C: THE COEFFICIENTS af ft CORRESPONDING TO RIGHT- 
HANDED 4- QUARK OPERATORS 

In (27), the expressions for effective parameters o^ 12 ^ corresponding to right-handed 
4-quark operators are 



«f (W 2 ) = 



Ci±i C F a s ( h 



N r 4tt 



(CI) 



with Ni(y 2 ) = 1 for i = 11, 12. For o^( 14 \, one should replace q by q, and have iVi3(V2) = 1, 



-^14(^2) = 0. V^(V 2 ) account for vertex corrections, H^{ViV 2 ) for hard spectator interactions 
with a hard gluon exchange between the emitted meson and the spectator quark of the B 
meson and Pi(V 2 ) for penguin contractions. The vertex corrections read 



V?(V 2 ) = { 



( f 1 




/ dx $y 2 


[x) 


Jo 








< / dx® V2 


(x) 


Jo 








f dx$ V2 


» 


, Jo 





(i = 11,12), 
(i = 13), 
(i = 14), 



(C2) 
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and 



( rl 



dx $6 2 (x 



V?(V 2 ) - < / dx$„Jx 



121n^ - 18 + g T (x)} , 
p J 

" _ 121ii — + 6-^r(l-x) 
p 



0, 



(i = 11,12), 
(i = 13) 
(< = 14), 



where $y(a;), $„(x), $ a (x), $fe(x), g(x), h(x) and gr(^) are defined in [31] and [15]. 
(V1V2) have the expressions: 



H^VM) = H° 12 (ViV 2 ) = 



fafviM f 1 J *f (p) 



(bvi,v 2 ) 







dp- 







/w 1 

Jo Jo 



X 



X (BVuV 2 )J p 



fdvf 
Jo Jo 



X 



du 



<f> Vl (v)$v 2 {u) rVi $ Vl (v)<!> V2 (u) 



uv 



uv 



HUV1V2) = and 

H+(y 1 v 2 ) = H+(y 1 v 2 ) 

H+(YiV 2 ) 

HUV1V2) 

where 



x (BVi,V 2 ) 



[ dv [ 

Jo Jo 



dp- 



du 



UW 



f_BfyJv2 f 

Vl ' V2) Jo 



P Jo Jo 



dv du 



UV' 



(C3) 



(C4) 



(C5) 



2X (BV 1: V 2 ) J Q 



dv / du 



P Jo Jo 



uuv 



-, (C6) 



X, 



(bVi,V 2 ) 







/v 2 
2m 



/ 2 2 



Xf™ = -fv.mv 2 



rn 



(C7) 



with q = Pb — PVi = Pv 2 - Here $f (p) is one of the two light-cone distribution amplitudes of 
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the B meson [48]. P^' p are strong penguin contractions. We obtain 



Pif{V 2 ) 



CfCXs 

4ttN c 

i 

i=u 

C F a s 
~ AnN c 

b 



C12 



^± i„!^ _ (n/ _ 2 )Gy 2 (0) - GV 2 (a c ) - G Va (l) - J 



5ln^-GV 2 (0)-Gy 2 (l) + ^ 



3 
(C8) 



(Cl2 



+ Cl 4 )^ 



(n / -2)G y2 (0)+Gv 2 (s c ) + G' y2 (l) 



GV a (0) + G Va (l) 



(C9) 



p^'P = p^'P = p+' p = P 1 4' p , where = mf/ml and the functions Gm 2 {s) and GW^s) are 
given by 



Gm 2 ( s ) — ~~ 4 / rf-u$y 2 (-u) / (irxa; ln(s — uxx — ie) 
Jo I Jo 



G v 2 { s ) = -4 / du$ V2 (u) 



dx xx ln(s — -uxx — ie) 
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